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ABSTRACT 


An analytical and experimental investigation has been made to 
evaluate a technique by which the performance of charring ablators in 
high enthalpy environments can be determined by testing in low enthalpy 
environments. Relationships are presented which indicate that proper 
simulation of a high enthalpy environment can be achieved by adjusting 
the heat input and stream oxygen content in the low enthalpy environment. 

Tests with a high and a low density phenolic nylon charring 
ablator were made at two enthalpy levels to validate the simulation 
technique. Values of the principal test parameters in the high enthalpy 
environment and the low enthalpy test environment were, respectively, 
enthalpy, 1000 and 5000 Btu/lbm, cold-wall convective heating rates, 

170 and 160 Btu/ft^ sec, stream oxygen mass fraction 0.23 and 0.05 • A 
comparison of the ablation performance data indicates that satisfactory 
simulation was obtained. 

Although the present tests were conducted at constant enthalpy 
heating, a method is outlined to extend the present technique to simulate 
the ablation performance at variable enthalpy in a constant enthalpy 
environment by step changes in oxygen concentration. 



CHAPTER I 


INTRODUCTION 

The charring plastics have been shown to provide a practical 
method for protecting the interior of a spacecraft from the atmospheric 
heating encountered during entry into planetary atmosphere. The behavior 
of this type of material when subjected to heating is very complex, and 
extensive experimental (Refs. 1 and 2) and analytical investigations 
(Refs. 5, and 6) have been undertaken to study this problem. 

Experimental studies have been conducted in a variety of test 
facilities which produce a wide range of environmental conditions. Many 
of these facilities provide the appropriate heating rate, but are charac- 
terized by low enthalpy levels. The high enthalpy facilities are usually 
restricted to tests of snail specimens and are therefore useful only for 
materials evaluation tests . Some of the low enthalpy facilities can be 
used to test much larger specimens and combined thermal and mechanical 
performance can be investigated. However, ablation performance at low 
enthalpy levels may differ from that which would be obtained at high 
enthalpy levels . 

In the present paper the conditions under which the performance 
of a charring ablator in a high enthalpy environment, typical of reentry, 
can be simulated in a low enthalpy test facility are examined. The 
results of Reference 7 are used to determine quantitative relationships 
between results obtained under different test conditions . Equations are 
presented to determine the test conditions which should be used in a 
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facility having a given enthalpy to produce the same material performances 
that would be obtained in an environment having a different enthalpy. 

Experimental results for a low-density and a high-density charring 
ablator are presented to determine if simulation of the ablation per- 
formance was achieved. 

The basic equations resulting from this study and a set of quali- 
fying test data are presented in Reference 8. The data of Reference 8, 
however, are preliminary and indicate only that the overall performance 
can be duplicated . In the present study the complete reproduction of 
the transient behavior is considered as well as the overall ablation 


performance. 



CHAPTER II 


APPROXIMATE ANALYSIS OF CHARRING ABLATORS 

An analytical model of the charring ablator is presented in 
Reference }. This model is shown schematically in Figure 1. The outer 
surface of the char layer is subjected to aerodynamic heating. As the 
material is heated, the material undergoes degradation, forming a char 
layer which provides both insulation and a high-temperature outer 
surface for reradiation. The heat passing through this layer is 
partially absorbed by pyrolysis at the interface between the char layer 
and the uncharred material, and the remaining heat is conducted into 
the uncharred material. The gases generated by pyrolysis transpire 
through the char layer and are injected into the boundary layer. The 
gases are heated as they pass through the char, which reduces the 
quantity of heat conducted to the pyrolysis interface. When these 
gases are injected into the boundary layer, the convective heat transfer 
is reduced. In addition to the gases produced by pyrolysis, a residue 
remains at the interface and adds to the thickness of the char layer. 

At the same time, char removal may occur as a result of thermal, 
chemical, or mechanical process. Thus the total char thickness may 
increase or decrease depending on the relative rates of formation and 
removal of the char. These processes are related quantitatively in the 
following section. 
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Figure 1.- Schematic diagram of charring ablator. 
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Governing Equations 

It is assumed that thermal properties of the char are functions 
only of temperature, that all heat flow is normal to the surface, and 
that gases transpiring through the char are at the same temperature as 
the char. The differential equation governing the temperature in the 
char is (from Ref. 3)> 


S /. St\ • — St ST /, n 

Sx\ Sx ) ^ Sx = pCp St (1) 

To solve this equation, conditions on the temperature or heat 
flux and mass transfer must be given at each boundary. In a charring 
ablator with low rates of mass transfer into the boundary layer and with 
char removal by oxidation only, the conditions at the char surface are 
as follows (Ref. 3), 

^c,net + 

Heat input Reradiation Conducted to 

interior 

where, for diffusion controlled oxidation with a Lewis number of 1, 


m c = AC e 


Qc,net 
(h e - hy) 


( 3 ) 


^c,net 



<V"c 



(M 


and 
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Pyrolysis occurs in a zone between the char and the uncharred 
material. In many cases this zone is very thin and can be treated as a 
plane. Reference 9 indicates that treating the pyrolysis zone as a 
plane is a reasonable assumption. In this case, the temperature of the 
char and the uncharred material must be the same at the pyrolysis plane 
and, in general, this temperature will be a function of the rate of 
pyrolysis, 

^char “ ^uncharred = F(mp) 

From an energy balance at the pyrolysis plane, 


- k 


St\ 

bxj 


char 


* iy%, 


- k 



uncharred 


(5) 


Solution of the above equations, together with the transient 
heat conduction equation in the uncharred material, provides an estimate 
for all aspects of the performance of a thermal protection system based 
on a charring ablator. Detailed application of these equations is 
presented in Reference 3. 


Approximate Equations 

An approximate analysis of the performance of charring ablators 
is presented in Reference 7 • The analysis is based on a number of 
simplifying assumptions . These assumptions are discussed in Appendix A. 
The analysis provides equations for calculating the amount of material 
degraded by heating at any time. The char thickness and the amount of 
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char removed can also be determined. However, the internal temperatures 
cannot be determined on the basis of such a simplified analysis. 

Qv* 

For the case when enthalpy, , and stream oxygen concentration 

Hco 

are constant and the oxidation process is diffusion controlled, the 
weight of char accumulated for the convective heat input, Qc, is (from 
Ref. 7), 



where 


a = C e A 


(7a) 


0 = 


1 + T|(aAh + h e - h^) 

Ah eff (l + aT l) 


(7b) 


^rP 

ae( + T^)(T r + T p ) 


(7c) 


and 


' Ah + ^e ~ *Hr\l - f + ^r (h e - hy) /l - f + arA ^ ^ 
Z^eff ^eff / f co ^eff V af / 


( 6 ) 
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The extremum char weight, W c m , is the weight of char which would 
accumulate if a specimen of infinite thickness was heated indefinitely 
at a given heating rate. 

The amount of material pyrolyzed is (Ref. 7), 



f(l + atj) 

and the total weight of material degraded is 

w < ■ + ( r 7.J0JH.) 

w + w. 

1 - f + aij 


(9) 


(10) 


(1 + aij) 


CHAPTER III 


SIMULATION OF ABLATION PERFORMANCE 

The equations presented in the previous section can be used to 
determine the performance of charring ablators exposed to various 
thermal environments with constant enthalpy. In this section, the 
conditions under which the same ablation performance can be obtained in 
different but constant enthalpy heating environments will be investi- 
gated. Relationships are determined that specify the required environ- 
ment for simulation. 

It is assumed that the enthalpy level of the test stream is 
relatively constant and much lower than the stream in which the per- 
formance is to be determined. This is the characteristic situation 
which is encountered when an attempt is made to evaluate reentry per- 
formance using large-scale specimens. 

An examination of the equations in the previous section shows 
that environmental parameters affect performance directly only at the 
heated surface. Therefore, any two environments so related that the 
inputs at the surface are the same should result in the same ablation 
performance of the material That is, if the surface has a given 
history of location and temperature, all phenomena beneath the surface 
are independent of the particular environment which produces the surface 
conditions . 
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The aerodynamic innuts to an ablating surface may be classified 
as thermal, chemical, and mechanical, and represented analytically as 
follows : 

Thermal 


"Jr + ^c,net ( lla ) 

Chemical 

m(0 2 ) (lib) 

Mechanical 

Pyj? + r^t (11c) 

Only oxidation is considered here under chemical inputs. In some cases, 
stream species other than oxygen may interact with the surface, bat the 
following procedure c^n be extended directly to such cases and they will 
not be considered here. 

No satisfactory general theory of the effects of mechanical 
inputs on the performance of charring ablators is available. However, 
experimental data indicate that at low shear stress levels, such stress 
has no effect on thermal performance and can be neglected. 

Both convective and radiative heating can exist in the environ- 
ments considered here. Either convection or radiation is the mode of 
heat transfer in most low enthalpy test facilities. The present 
discussion will be limited to the case where convective heating only 
is present in the low enthalpy environment, and both convection and 
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radiation may exist at the high enthalpy environment. Appendix B 
discusses the case when radiation is the prime mode of heat transfer in 
the low enthalpy environment . 

Two conditions that should exist to obtain simulation are: 


(W c ' 1 = (W c ) 2 and (W c>m ) 1 = (W c , m ) 2 (12) 

The ablation performance in environment 2 is to be simulated in environ- 
ment 1. Assuming convective heating to be the only mode of heat 
transfer in environment 1 and equating the extremum char weights, 
equation (8), for the two environments , 

f Ah + h e - h v \l - f 

£h eff a£h eff J f 

( 13 ) 




Since the material properties in equation (13) are considered invariant 
with respect to the environment, equation (13) simplifies to 


C e 'l 

1 / C e \ 1 

\h e “ hw; 

1 

rt> 

1 

r 

ro 

«flr (1 - f + an) 

1 + ^=o 1 - * 


CM 


The material and the environment interact directly only at the 
interface between them. Therefore , the net aerodynamic heat input 
should be the same in the two conditions, that is. 
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^c,net^ “ ^c,net + ^r^ 

To obtain the heat input needed for simulation, equation (6) is first 
used to calculate the char weight V c ^ produced in environment 2 for 
a given convective heat input. This char weight, W C2 , which equals 
W Cl , and the known environmental parameters of the first environment 
are then used in equation (6) to determine the convective heat input 
required in environment 1 to obtain simulation. 

Because of the transcendental character of equation (C), 
equations (6) and (10) cannot be solved explicitly to obtain a direct 
relation between weight of material degraded and heat input. Therefore, 
it was found advantageous to solve these equations for the required 
simulation test conditions, as well as char weight and weight of 
material degraded, on a digital computer. 

To summarize, the following procedure may be used to obtain the 
test conditions required to simulate, in a known constant enthalpy 
ground test environment, the ablation performance of a charring ablation 
in a different but constant enthalpy environment: 

1. Equation (l4) relates the oxygen concentrations of the 
two environments each of different but known enthalpies. 

2. The convective heat input required for simulation is obtained 
by using equation (6) to determine the char weight produced by the 
second environment $ then by utilizing the identity in equation (12) and 
all of the known environmental parameters, equation (6) is solved for 
the heat input for the ground test. 
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The weight of char formed, equation (6), and the weight of 
material degraded, equation (10), depend on an integrated heat input 
and are independent of the heating rate history. Therefore, for tests 
to evaluate degradation characteristics, test time is related to flight 
time, through equation (15), 


t 


1 



+ aq r ) 2 dt 


^c,net^ 


(16) 


To reproduce transient effects, however, test times must be the same. 


CHAPTER IV 


TEST SPECIMENS 
Materials 

Two char forming ablation materials were used in this investiga- 
tion. One material contained, by weight, 50-percent powered nylon and 
50-percent powered phenolic resin and was pressure molded to a density 
of 75 lb/ft^. The other material contained, by weight, 4o -percent 
powered nylon, 57 -percent powered phenolic resin, and 25-percent 
phenolic microspheres, and was pressure molded to a density of 5^ lb/ft^. 
These materials represent two types of ablation materials that have been 
extensively tested. The low-density material is representative of a 
group of materials that have been flight tested and have application to 
future spacecraft. 

The thermophysical properties for these materials, required by 
the present analysis, are given in Table I. These properties were 
obtained from Reference 10 and unpublished data. The thermal transport 
parameter, 0, as defined by equation (7c), was evaluated for the present 
study using a reference temperature equal to the radiation equilibrium 
temperature based on the nominal test conditions, that is. 



Material property 
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Ah, Btu/lbm 
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- MATERIAL PROPERTIES 


Low-density 
phenolic nylon 

High -density 
phenolic nylon 

0.115 

0.19 

5,000 

5,000 

0.75 

0.75 

0.53 

0.6 

4,3^0 

4 , 3^0 

0.6 

0.6 

0.9 

0.9 

3.61 x 10" 4 

3.19 x 10" 4 

16 

30 
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The reference conductivity for each material was obtained with the 
relationship, 



where k is the temperature dependent char conductivity for the 
respective material from Reference 10. 

Geometry 

The test specimens consisted of 2 .5 -inch -diameter cylinder, with 
a 4 .7-inch radius spherical face, 1.25 inches thick, machined from large 
blanks of the phenolic nylon materials. The specimen geometry provides 
uniform heating and pressure across the front surface (Figs. 2 and 5)* 

These cylinders were bonded to bakelite holders (Fig. 2). There 
were eight specimens of each material, two of which contained 1-inch- 
diameter instrumented plugs. 

Bach plug contained four No. 3 6 chromel-alumel thermocouples 
(Fig. 3). The thermocouple junctions were located nominally at 0.124, 
0.25, 0.375, and 0.5 inch from the front surface. The thermocouple 
installation recommended in Reference 11 was used. The actual thermo- 
couple bead locations were determined by X-ray and are noted on the 
corresponding temperature histories. 
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Instrumented Specimen 


(a) Low-density phenolic nylon 


Figure 2.- Typical specimens before testing 
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(b) High-density phenolic nylon. 


Figure 2.- Concluded 
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Figure 3.- Sch emat ic of typical instrumented specimen. 



CHAPTER V 


CONDITIONS AND PROCEDURES 
Test Facilities 

The 2 -inch supersonic arc tunnel at the Langley Research Center 
was used for the tests conducted at the low enthalpy conditions. This 
facility is described in Reference 12. A conical nozzle with a throat 
diameter of 1.5 inches and an exit diameter of 4 inches was used. The 
desired stream oxygen concentrations are obtained by mixing air and 
nitrogen . 

The high enthalpy direct current plasma generator at the Langley 
Research Center was used for the tests conducted at the high enthalpy 
conditions. This facility is described for a subsonic configuration 
in Reference 12. For the present study the facility was operated in 
a supersonic configuration with a 1-inch-diameter minimum and a 4-inch- 
diameter exit. Stream oxygen concentrations were obtained by mixing 
oxygen and nitrogen. 


Test Conditions 

The nominal test conditions for the environments were as follows 


Total enthalpy, Btu/lbm 

5,000 

10,000 

Free-stream Mach number 

3 

3 

Model stagnation pressure, atm 

0.37 

0.04 

Stagnation heat-transfer rate, Btu/ft^-sec 

140 

155 

Free-stream oxygen content, percent 

6 

23 

Model test time, sec 

varied 

varied 
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These test conditions were determined in the following manner. 

The high enthalpy test condition was established first by calibration 
test;* using the pressure and heating rate probes of the same size and 
shape as the specimen. Then assuming an enthalpy of 3 >000 Btu/lbm, the 
corresponding heat input and stream chemical composition required for 
simulation were determined using equations (6) and (14-) and the inputs 
given in Table I. Convective heating was the only mode of heat transfer 
present in these two environments. Tne hot-wall enthalpy used in these 
calculations was based on the radiation equilibrium temperature for the 
test condition being considered, that is, 



Models were tested for 30, 6o, 90, and 120 seconds. Tne instru- 
mented models used to determine transients simulation were tested for 
120 seconds. Therefore, each model tested at the high enthalpy condi- 
tion for a given time had a corresponding model tested at the prescribed 
low enthalpy condition for the same length of time. 

The present analysis is restricted to char mass loss fc„. diffusion 
controlled oxidation only. Therefore, it was necessary to sstablish 
the oxidation regime the present tests would provide. The detailed 
ablation analysis of Reference 3 was used to compute the ratio of the 
char mass loss rate to the uiffusion controlled mass loss rate as a 
function of cold-wall convective heat rate. The results of these 
calculations are shown in Figure **•. As can be seen, the present tests, 
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with a nominal cold-wall heating rate of 150 Btu/ft^-sec, will provide 
a diffusion controlled oxidation regime. 

Pressure distribution .- Stagnation pressures and pressure distri- 
butions were made on a pressure model (Fig. 5(a)) having the same 
external shape as the test models (Fig. 3)* The pressure distributicns 
over the model in both facilities are shown in Figure 5(b). The same 
distribution was obtained in both test streams . Stagnation-pressure 
measurements were obtained before and after each model was tested. 

These pressures were measured with an accuracy of ±10 percent. 

Heating rate .- Stagnation heating rates and heating rate distri- 
butions were measured with an 0.032-inch-thick inconel calorimeter 
model (Fig. 6(a)) having the same exterior shape and dimensions as the 
test models (Fig. 3)* The heating rate distributions over the model in 
both facilities were essentially the same (Fig. 6(b)). Figure 6(b) 
shows that the heating rate over the 1 -inch-diameter instrumented plug 
was essentially constant . This would provide one-dimensional heating 
in the region on the test models where measurements were made. The 
measured heating rates are believed to be within ±10 percent of the 
actual heating rates. 

Enthalpy . - Total enthalpy of the test streams was determined 
using References 13 and 14 and the measured stagnation pressure and 
heating rate for each model. The enthalpies are given in Table II for 
each model. Hie accuracy of these values is ±10 to 20 percent due to 
the dependence of the value on the heating rate and pressure. 


SSURE ORIFICES 



Figure 5*- Pressure probe end pressure distribution. 









O High Enthalpy 
□ Low Enthalpy 



(b) Heating rate distribution. 
Figure 6.- Concluded. 
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Test Procedure 

For each test, the facility operating conditions determined from 
previous calibration tests were established. Heating rate and pressure 
measurements were made prior to and after the model was exposed to the 
test stream for a predetermined time. Periodic checks were made on the 
test conditions and adjustments were made as required. All temperature 
and pressure data were recorded by the Langley control digital data 


recording facility. 



CHAPTER VI 


RESULTS AND DISCUSSION 

Eight pairs of specimens were tested in lov and high enthalpy 
environments related through equations (6) and (l4). Both specimens of 
a pair were tested for the same length of time. Four of the pairs of 
specimens were made of a low-density phenolic nylon material, each 
specimen has a LD prefix to the specimen number, and four pairs were 
made of a high-density phenolic nylon material, each specimen has an 
HD prefix to the specimen number. The test data are tabulated in 
Table II . All thickness measurements were made with an accuracy of 
±0.01 inch. 


Low-Density Phenolic Nylon Specimens 
Photographs of the tested specimens are shown in Figure 7- For 
each specimen tested at high enthalpy there is a corresponding specimen 
tested at low enthalpy. The total recession, char thickness, and thick- 
ness of material degraded for each specimen are given in Table II and the 
data are plotted in Figure 8. Also shown in Table II are the required 
test conditions at low enthalpy specified by the present analysis and 
the actual test conditions. As can be seen the actual test conditions 
were very close to the required test conditions. 

From a visual comparison of the paired specimens in Figure 7, it 
is seen that the char surface, thickness, and char structure are similar. 
The cracked pattern visible on the surface of each specimen is believed 
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TABLE II.- SUMMARY OF TEST DATA 



(a) Instrumented specimens 
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Figure 7 *- Continued. 
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(f) Specimen LD-5 tested 90 seconds at low enthalpy. 
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to be formed during specimen cooling. Motion pictures taken during 
testing do not show this pattern in either environment. 

The surface recession, char thickness, and material degraded 
thickness data from Table II are plotted in Figure 8 as a function of 
time. The agreement between these thicknesses for each pair of speci- 
mens is considered satisfactory since the difference is within the 
accuracy of the measurements, ±0.01 inch. Figure 8, therefore, indicates 
that the same ablation performance was obtained in the two related 
environments . 

Further indication of the simulation of the ablation performance 
can be obtained from a comparison of the internal temperature histories. 
Figure 9 is a comparison between the internal temperatures of a model 
tested in the high enthalpy enviroranent and the low enthalpy environment. 
As can be seer '-he agreement is very good . 

An examination of the data presented in Figures 7, 8, and 9 shows 
that the transient ablation performance in a high constant enthalpy 
environment of the low-density charring ablator considered here was 
satisfactorily reproduced in a low constant enthalpy test stream speci- 
fied by the present analysis. 

A comparison between calculated and measured char thickness and 
material degraded thickness for each of the low-density material 
specimens is shown in Table III. The results predicted by the present 
analysis are in poor agreement with the test results. The agreement 
between measured and calculated material degraded thickness is not as 
bad as that between measured and calculated char thicknesses. 
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Figure 9*- Internal temperature histories for lev -tensity phenolic nylon specimens 
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Specimen 

No. 


LD-15 

l4 


13 

7 


l6 

5 


3 


HD-10 

9 


8 

7 


13 

11 


3 

6 


Test 

data 

Present analysis 

Char 

thickness, 

in. 

Material 

degraded, 

in. 

Char 

thickness, 
in . 

Material 
degraded, 
in . 

0.15 

0.17 

0.172 

0.204 

0.14 

- 

1 

0.17 

0 1 72 

0.204 

0.2 

0.31 

0.278 

0.344 

0.2 

0.31 

0.278 

0.344 

0.22 

0.4i 

0.315 

0.480 

0.23 

0.42 

0.315 

0.480 

0.23 

0.52 

0.420 

0.546 

0.21 

0.49 

0.420 

0.546 

0.07 

0.08 

0.078 

0.099 

0.08 

0.10 

0.078 

0.099 

1 - ■ 

0.12 

0.17 

0.129 

0.173 

0.12 

0.18 

0.129 

0.173 

0.16 

0.23 

0.165 

0.206 

0.17 

0.25 

0.165 

0.206 

0.18 

0.33 

0.196 

0.294 

0.2 

0.35 

0.196 

0.294 


Reference 15 


in. 


0.21 


in. 


0.18 

0.18 


0.275 

0.29 


0.375 

0.385 


0.48i 

0.494 
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These differences may be accounted for in the following manner. In the 
present analysis it is assumed that the heating rate is constant over 
the time period. Figure 7 shows that the specimen geometry approaches 
a hemisphere, with time, resulting in an increase in heating with time. 
This increase would cause a decrease in char thickness due to an 
increase in surface recession. 

Another phenomenon present in these tests that was not accounted 
for in the analysis is the lateral flow of internal gases. This effect 
is discussed in detail in Reference 15. Also, Reference 15 presents 
test data, for a low-density phenolic nylon material similar to the 
material in the present study, that were successfully matched only when 
the lateral flow of internal gas was considered. Briefly, this lateral 
flow results in a reduction in mass transfer into the boundary layer 
caused by the lateral flow of pyrolysis gases and the inflow of gases 
at the surface. This lateral flow results from a combination of c u .r 
permeability, external pressure field, and the ratio of char thJc*uiess 
tc specimen body radius. Incorporating the lateral flow analysis of 
Reference 15 into the transient ablation analysis of Reference 5 and 
applying this to the present tests gives the indicated results shown in 
Table III. In addition, it was found necessary to reduce the char 
density from 16 lb/ft^ to 12 lb/ft^ to obtain the results in Table III. 
A measured char density of 16 lb/ft^ was used to determine the nominal 
test conditions, however, a char density of 12 lb/ft^, which was also 
measured, is reported In Reference 10 for this material. These chars 
were formed under different conditions. As can be seen, the agreement 


^0 


between calculated and measured char thickness and thickness of material 
degraded is satisfactory. Calculations indicated that the internal flow 
phenomenon effected the test results by about the same magnitude in both 
high and low enthalpy environments, see Table III. 

High -Density Phenolic Nylon Specimens 

Photographs of the tested specimens are shown in Figure 10. For 
each specimen tested at high enthalpy there is a corresponding specimen 
tested at low enthalpy. The total recession, char thickness, and thick- 
ness of material degraded for each specimen are given in Table II and 
the data are plotted in Figure 11. Also shown in Table II are the 
required and actual test conditions for the low enthalpy tests . 

From a visual comparison of the paired specimens in Figure 10, 
it is seen that the char surface, thickness, and structure are similar. 
As was previously mentioned, the cracked pattern present on all the char 
surface is believed to occur during cooling. The source of the sunburst 
pattern on the specimen tested for 30 seconds at high enthalpy is not 
known. 

Figure 11 shows the pyrolysis zone and outer surface location 
histories for the high density phenolic nylon specimens . Again the 
differences are within the accuracy of the measurement i . The agreement 
is very good, indicating satisfactory reproduction of the transient 
surface and pyrolysis zone behavior. 

Figure 12 shows internal temperature histories for models tested 
at the low and high enthalpy conditions . The agreement is satisfactory. 




INCHES 





Figure 10.- Continued. 




44 



(h) Specimen HD-6 tested 120 seconds at low enthalpy. 





Figure 12.- Internal temperature histories for high-density phenolic nylon 


The differences may be attributed to the heating rate experienced by 
the specimen tested at low enthalpy which was higher than that required 
for simulation. 

An examination of the data presented in Figures 1C, 11, and 12 
shows that the transient ablation performance in a high enthalpy 
environment of the high-density charring ablator considered here was 
satisfactorily reproduced in a low enthalpy test stream specified by 
the present analysis. 

A comparison between measured char thickness and material degraded 
thickness and the corresponding thicknesses computed by the present 
analysis for each specimen is shown in Table III. The agreement is 
satisfactory. Hence, it is assumed that the lateral flow of internal 
gases present in the low-density phenolic nylon specimens was not present 
in this material for these tests. Under some test conditions, however, 
the lateral flow of internal gases may be important for this material. 

In summary, relationships have been obtained that determine the 
heat input (eq. (6)) and the stream oxygen content (eq. (lh- ) ) required 
in a low constant enthalpy environment to simulate the ablation per- 
formance of a charring ablator in a high constant enthalpy environment. 
Two types of phenolic nylon ablation materials have been tested in 
environments related by equations (6) and (l4) . A comparison of these 
test data (the char recession, char thickness, thickness of material 
degraded, internal temperatures histories, and physical appearance) will 
indicate if satisfactory simulation of the ablation performance can be 


obtained . 



A brief comparison of the test data is given in Table IV. Hie 
reproducibility at low enthalpy of the char recession, char thickness, 
and thickness of material degraded measured at high enthalpy, is deter- 
mined by the percentage difference between these thicknesses at low and 
high enthalpy . The ratio of the accuracy of the measurement, ±0.01 inch, 
to the thickness measured is also given in Table IV as the percentage 
of accuracy. All percentages are based on the thiv .messes measured at 
high enthalpy. It is believed that any significant difference in 
ablation performance could be discernible by the accuracy used for this 
study. Table IV, therefore, indicates that the thickness data from 
these two sets of tests are in good agreement. 

Hence, a comparison of the physical appearance and the internal 
temperature histories, together with the thickness data, indicate that 
satisfactory simulation of the ablation performance of these two 
materials was obtained. 


Variable Enthalpy Heating 

The present analysis applies to the simulation of the performance 
of a char forming ablator subjected to high constant enthalpy heating in 
a low constant enthalpy heating environment. The preceding test results 
indicate that this analysis can be satisfactorily used to obtain good 
simulation. In this section a method will be outlined that may be used 
to simulate the performance of a material subjected to high, variable 
enthalpy heating, typical of entry into a planetary atmosphere, in a 
constant enthalpy heating environment. 
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TABLE IV.- COMPARISON OF THICKNESS DATA 


Time, sec 

Low-density material 

High -density material 

30 

6o 

90 

120 

30 

6o 

90 

120 

Recession, 









Percent difference 

25 

0 

0 

3 

100 

20 

14 

7 

Percent accuracy 

25 

8 

5 

3 

100 

20 

14 

7 

Char thickness. 




■ 





Percent difference 

8 

0 

5 

H 

14 

0 

6 

5 

Percent accuracy 

8 

5 

5 

■ 

l4 

8 

6 

5 

Material degraded. 







■ 


Percent difference 

0 

0 

2 

6 

24 

6 

D 

6 

Percent accuracy 

6 

3 

2 

2 

12 

6 

D 

3 
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Figure 15 shows a typical enthalpy and cold -wall stagnation 
heating rate history for a manned lifting-body vehicle entering the 
earth's atmosphere. To apply the approximate ablation analysis of 
Reference 7, a constant enthalpy heating condition must exist. 

Therefore, an approximation of the enthalpy history in Figure 13 is 
needed. To approximate constant enthalpy heating, a variable enthalpy 
heating trajectory, such as the one in Figure 13, may be divided into a 
number of steps or intervals over which the enthalpy may be considered 
constant. A typical step is indicated in Figure 13 between 500 and 
600 seconds. It may, therefore, require a large number of steps to 
approximate the enthalpy history. 

Applying the approximate ablation analysis co the approximated 
enthalpy history and undisturbed heating history, the char weight 
formed during each step is obtained from equations (6), (7), and (8). 

Hie total weight of material degraded during a step is obtained from 
equation (10) . The heat input for a given step is obtained by inte- 
grating the heating rate over that step (Fig. 13) • The initial weights 
used in equations (6) and (10) are the final weights from the previous 
step. 

Hiis approximation to the enthalpy history allows the present 
analysis to be used to simulate the performance of a charring ablator 
subjected to a varying enthalpy heating environment in a. constant 
enthalpy environment . The required oxygen concentration of the test 
stream for simulation during a given step is obtained from equation (1^ ) . 
Since the char veight and the extremum char weight are the same in flight 



led lifting-body entry vehicle. 


§2 


and test for each step, the magnitudes of W c and W c>m are obtained 
from equations (6), (7), and (8) using the approximated flight condi- 
tions. The test heat input is found by scaling the flight heat input 
through equation (6) . The test time for a step may be obtained from 
equation (l6). Since the values of C e , Qq, and t obtained with 
equations (6), (14) , and (l6) apply to only one step, the above process 
must be repeated for each step using the final weights from the previous 
step as the initial weights for the next step, therefore, the test 
required to simulate the ablation performance at varying enthalpy will 
consist of a series of steps of different ^ e , , and t at constant 


enthalpy . 



CHAPTER VII 


CONCLUDING REMARKS 

The conditions under which the performance of a charring ablator 
in a high enthalpy environment can be simulated in a low enthalpy 
environment have been examined. An epproximate ablation analysis has 
been used to obtain quantitative relationships between a low enthalpy 
heating environment and a high enthalpy heating environment to obtain 
the same ablation performance of a charring ablator in both environments . 

Tests were made to determine the ablation performance of a low- 
density phenolic nylon material and a high -density phenolic material in 
air at a constant high enthalpy heating condition. Identical specimens 
of each material were also tested in a constant low enthalpy test stream 
related to the high enthalpy environment by the aforementioned relation- 
ships. These tests indicated that satisfactory simulation of the 
ablation performance of these two materials was obtained. 

The low-density phenolic nylon material performance was affected 
by its high char permeability and external pressure environment. 

However, the material performance appeared to be affected to the same 
degree in both high and low enthalpy environments. The high -density 
phenolic nylon material did not appear to be affected in this way. 

The testB presented in this paper indicate that satisfactory 
simulation was obtained. However, convective heating at constant 
enthalpy was the only case considered. A method is outlined to extend 
the present technique to simulate the ablation performance at variable 
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enthalpy with convective and radiative heating in a constant enthalpy 
environment with convection only. Simulation of the latter case at 
constant enthalpy with radiation as the prime source of heat is also 


discussed. 
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APPENDIX A 


ASSUMPTIONS 

A number of simplifying assumptions must be made if approximate 
analytical relations governing the performance of charring ablators 
are to be obtained. Hie assumptions used in the analysis on which the 
present paper is based sure: 

1. Hie specific heat of the char is neglected. This assumption 
reduces the problem from a transient to a quasi-steady-state condition. 
Calculations based on this assumption are compared to numerical solutions 
which include specific heat effects in Reference 16 and the agreement is 
satisfactory. 

2. The thermal conductivity of the char is assumed to be a 
function of temperature as follows: 

^r(^j + lg)(T s *p) 

(Tp + T§)(T r + T p ) 

The thermal conductivity of a porous material such as a char layer will 
increase with increasing temperature because of radiation across the 
void spaces. Relatively large changes In heating rate are required to 
obtain significant changes in the surface temperature because the 
surface reradiation energy at a rate proportional to the fourth power 
of the temperature. Therefore, by proper selection of the reference 
temperature, T r , the error introduced can be minimized. 
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3 . Char removal is assumed to occur only as a result of a 
diffusion limited oxidation mechanism. This yields a maximum rate at 
which char can be removed by oxidation, and will be conservative if 
oxidation is the only mechanism of char removal . 

The product of combustion is assumed to be carbon monoxide, 
and all of the heat generated by the reaction is assumed to go into the 
surface. 

5. For moderate mass transfer the reduction in convective 
heating resulting from mass transfer into the boundary layer is assumed 
to be directly proportional to the rate of mass transfer (Ref. 3) • It 
is also assumed that the char which is oxidized and the pyrolysis 
products which are injected into the boundary layer sure equally 
effective in blocking the convective heat input. 

6. The Lewis number is assumed to be unity, that is. 


«Le = 


Npr _ CppD 
»Sc = k 


1 


This indicates that the mechanism of oxygen diffusion is the same as 
heat transfer in the boundary layer. Although the value of the Lewis 
number is a function of the temperature and the gas being considered, 
the Levis number is roughly unity for most gases of interest. Also, 
Reference 17 indicates that heat transfer at the stagnation point is 
not sensitive to the effect of dissociation, that is, N^g = 1. 


APPENDIX B 


SIMULATION WITH RADIATIVE HEATING 

When radiant heating is the prime heating mode, it is assumed 
that the gas flow over the surface is used only to provide the required 
oxidation of the surface, and to remove the gaseous products of pyrolysis. 
The convective heating is approximately equal to zero, as is the 
difference between the enthalpy of the stream external to the boundary 
layer and at the surface. However, the ratio of the heating rate to 
the enthalpy difference has a definite value which is the heat transfer 
coefficient. 


9c, net 
“ h w 


h e 


(B-l) 


The following relations sure obtained from Reference 7, using 
equation (B-l): 


dW c 

dt 



ip - ay 


The heat balance at the surface becomes 


where 


q l + ^r 



(B-2) 


(B-5) 


q-L = ayAh 


6 0 
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Combining equations (B-2) and ( B-3 ) , 


owj r a>Ah 
+ 


Ah e ff ^eff 


<-> 


Hie extremum char weight then becomes, 


«c,m = P 


^ 1 - fV °^r + Ah \ 1 

\ f A a7Ah eff ^eff/ " 


(B-5) 


For the test environment, equations (6), (7b), and (10) become: 



where 


* = 


1 + 


r|aAh 


Ah e ff(l + an) 


(B-7) 


and 


W = 


W„ 


- W 1 + (lTaT[)^ 
1 - f + arj 


(B-8) 


1 + ar^ 
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Assume that both convective and radiative heating are present in 
the flight environment. Hie extremum char weight for the flight environ- 
ment is given by equation (8) . The oxygen concentration in the test 
stream is determined by equating equations (8) and ( B-5 ) , 

(^1 

(^e) 2 

2 



The required radiant heat rate is obtained from the equation 

far> x = far > 2 + fao^ 1 ‘ " ^c“c + V^P^] 

(B-10) 

The required test time is related, as before, to the flight 
time by equation (16). 

In summary, to simulate the ablation performance of a charring 
ablator in a flight environment, in which both convective and radiative 
heating are present, the required oxygen concentration in a test stream 
is determined from equation (B-9), the radiant heating rates determined 
by equation (B-10), and the test time by equation (16) . The test 
stream heat transfer coefficient required for simulation is obtained 
from equation (B-6) using the flight char weight from equation (6). 


